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Abstract 

Q>^' ' We show how to derive the complete set of electromagnetic corrections to 

0^ . the Nambu-Jona-Lasinio (NJL) model of the nucleon. Our results enable an 

accurate estimate of the electromagnetic contribution to the neutron-proton 
mass difference within this model. At the same time, our procedure demon- 
strates the way to calculate the complete set of meson corrections to the NJL 
model that maintains chiral symmetry. 

(N 
> : I. INTRODUCTION 

O I We have recently shown how to calculate all possible electromagnetic corrections, of order 

Q ■ e^, to any quark or hadronic model whose strong interactions are described nonperturbatively 

Qs ■ by integral equations |1[. Here we would like to apply our method to the three-quark NJL 

model of the nucleon and in this way derive the expression for the neutron-proton mass 

difference that is due to a complete set of electromagnetic interaction within this model. 

For this purpose it is useful to summarise the main model-independent results obtained 

^ ■ in Ref . |I[] . Within the framework of relativistic quantum field theory, the strong interaction 

C ■ Green function G describing a system of quarks or hadrons is given nonperturbatively by 

the integral equation whose symbolic form is 

^'- G = Go + GoKG. (1) 

c^ \ 

The complete set of lowest order electromagnetic corrections to the Green function G, de- 
noted by (5G, then follows from Eq. ([I|) on topological grounds: 

5G = 6Go + 6G0KG + GoSKG + G0K6G + {G^R-'G + G^KG" + GqK^'G'') D^^ (2) 

where SGq is the complete set of electromagnetic corrections to the free Green function Gq, 
and D^i, is the photon propagator that connects the appropriate currents (quantities with a /x 
or u superscript). Unlike 6G which has internal photons coupled everywhere, 6K consists of 
the strong interaction potential K with all possible photon insertions except those that start 
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or finish on an external quark or hadron leg. All currents G*^, Gq, and K'^ are constructed by 
the gauging of equations method 0] that effectively attaches external photons in all possible 
ways to the corresponding strong interaction quantities G, Gq, and K, respectively. Using 
this method one obtains 

G^' = GT^'G, r'^ = r(; + ir^, t^ = Go'g^Go'. (3) 

With the current G^ specified in this way, Eq. (Q) can be formally solved to give 

5G = GAG, A = 6K + Gq^SGoG^^ + (F'^GT'^ - T^GoT'^q) D^^. (4) 

The quantity A as given by Eq. (H) is the key result that describes the complete set of 
electromagnetic corrections to any observable of the strong interaction model in question. 
For example, if the strong interactions admit a bound state of mass M and wave function 
■0, then the complete set of lowest order electromagnetic corrections to M is given by 5M = 
il)Aip/2M. It is a feature of our approach that the gauge invariance of such electromagnetic 
corrections is a result of their completeness. 

II. NJL MODEL FOR THE NUCLEON 

The simplest NJL Lagrangian density C is defined in terms of the iso-doublet (two 
fiavours) colour-triplet quark field \1/ as 

C = ^{i^-mQ)^ + G [{^^f - (%5r^)2] (5) 

where r is the vector of isospin Pauli matrices, and ttiq is the bare quark mass [for itiq = 
Eq. (^) is chiral invariant]. The model of the nucleon considered here is described by the 
three-quark wave function that satisfies a four-dimensional Bethe-Salpeter (BS)-like three- 
body integral equation with pair interaction kernels given by the lowest order qq irreducible 
diagrams corresponding to the Lagrangian of Eq. (|^), namely 
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iG [{IsIfIc)^ iUfQj - (75T/c). ■ (75T/c),] (6) 



where Is, If and Ic are the unity operators in the Dirac, fiavour and colour spaces respectively, 
with the subscript i (j) indicating that the corresponding operators act in the i-th (j-th) 
quark's one-particle space. In Eq. @ and in the equations below we treat the quarks as 
distinguishable particles as the inclusion of antisymmetrisation can always be taken care of 
at the end. In the mentioned BS-like integral equation, the quark propagator d{p) satisfies 
the (nonlinear) Dyson- Schwinger equation 

d{p) = doip) + do{p)np)d{p) (7) 

where do{p) is the bare quark propagator and the dressing term S is taken in the so-called 
Hartree approximation: 

^^^^ = ^^/ (0ii^'^(^)^M - A'tT[d{k)A,]}. (8) 
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FIG. 1. The Bethe-Salpeter equation for the quark-diquark Green function. 

Here A^ is the Lorentz four- vector (Iglfl^'y^Tlc), and the trace "tr" is over the Dirac, 
flavour and colour indices. The kernel of Eq. (^) is effectively separable so that the original 
three-body equation for the three-quark system can be reduced down to a quark-diquark two- 
body equation, and it is this latter form which we shall use to calculate the electromagnetic 
corrections. In the channel where the two interacting quarks form a scalar, isoscalar, positive 
parity di quark with colour 3, the kernel of Eq. (H) reduces to 

where iii2 (/1/2) are triples of initial (final) quantum numbers of the first and second particle 
0. Then the diquark propagator is 

D{p) = ^ _ 2ig}i^p2) ^^^^^ n(p2)5,, = -/^tr[75r.d(p + A;)75r,d(A;)] (10) 

and the quark-diquark interaction kernel is given by the quark exchange term 

K{p',p)=^,d{p' + ph,f3'^f3^'. (11) 

With Go = dD in Eq. (|I|), the resulting equation for the quark-diquark Green function has 
the diagrammatic form illustrated in Fig. ||. 

III. ELECTROMAGNETIC CORRECTIONS TO THE NJL MODEL 

All electromagnetic corrections can be found by applying the general formulation of 
Eqs. (|l])-(|D to the particular case under consideration. Thus in the case of a single quark, 
even though the solution of Eqs. (0) and (|]) is known to be d{p) = i{^ — m)~\ where 
m is a constituent quark mass (which is not zero even if rriQ = 0), these equations are 
needed for the proper construction of the external and internal (with respect to the quark 
propagator) photon currents. By identifying Eq. (|^) as a special non-linear case of Eq. (|l]) 
(non-linear because S contains d) and gauging both Eq. (|^ and Eq. (^, we obtain that the 
electromagnetic current d^ of the dressed quark propagator satisfies the equation 



rf'^(p) = d{p + q) 



d^k 



r + iG j ——{K-d^^{k)K - A"tr[rf^(A;)A„]} 



d{p) (12) 



which is linear in d^ and can be easily solved [in Eq. (0) the momentum of the incoming 
photon, g, is contained implicitly in all the d^ functions]. The electromagnetic corrections 
to the dressed quark propagator, Sd{p), can then be found in a similar manner. 
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FIG. 2. The electromagnetic corrections to the NJL model. Shown are the diagrams making 
up the A of Eq. (|): (a) dK (first diagram) and Go^5GoGq^ (the rest), (b) diagrams of r^GF'^D^^ 
(also contributing are the last three diagrams with initial and final states interchanged), and (c) 
diagrams making up the subtraction term — rpGorQ^^jy. For lack of space we have not shown the 
corrections to the quark propagators in (a) due to diagrams such as (d). 

We can similarly identify Eq. (|T0|) for the diquark propagator D with Eq. (0) and therefore 
write down the diquark current D^^ and the electromagnetic corrections to the diquark 
propagator 6D as 

D^ = DWD; 5D = D6nD + DWDWDD^^. (13) 

Finally, we can write down the complete set of electromagnetic corrections 5G corre- 
sponding to the three-quark NJL model by identifying the quark-diquark equation of Fig. |l| 
with Eq. (|I|). The resulting expression for A is shown diagrammatically in Fig. 0. 

Clearly, our procedure for finding a complete set of electromagnetic corrections in the 
NJL model is of a general nature, and can be used to include complete sets of lowest order 
corrections due to other particle exchanges. Indeed we have recently applied our method to 
derive the complete set of pionic corrections to the three-quark NJL model (where previously 
only a part of such corrections were included 0]), as well as to some constituent quark models 
with confinement |§ (thereby clarifying some recent discussions regarding this matter IPJ^I)- 
As our pionic corrections are complete, axial current is conserved exactly in the limit of 
massless bare quarks. This feature is important for maintaining chiral symmetry in the next 
to leading order approximations. 
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